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We consider the behavior of massive Dirac fields on the background of a charged de-Sitter black 
hole. All black hole geometries are taken into account, including the Reissner-Nordstrom-de-Sitter 
^ ' one, the Nariai case and the ultracold case. Our focus is at first on the existence of bound quantum 

mechanical states for the Dirac Hamiltonian on the given backgrounds. In this respect, we show 
that in all cases no bound state is allowed, which amounts also to the non-existence of normalizable 
time-periodic solutions of the Dirac equation. This quantum result is in contrast to classical physics, 
and it is shown to hold true even for extremal cases. Furthermore, we shift our attention on the very 
interesting problem of the quantum discharge of the black holes. Following Damour-Deruelle-Ruffini 
approach, we show that the existence of level-crossing between positive and negative continuous 
energy states is a signal of the quantum instability leading to the discharge of the black hole, and 
in the cases of the Nariai geometry and of the ultracold geometries we also calculate in WKB 
approximation the transmission coefficient related to the discharge process. 
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I. INTRODUCTION 
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Black holes in de Sitter space are an interesting subject of investigation, both on the theoretical side and on 
, the experimental one. On one hand, the contextual presence of a black hole event horizon and of a cosmological 
event horizon, to be associated with the corresponding quantum emission of thermal radiation [l| is a feature which 
CNj | enriches the framework of black hole thermodynamics in itself and also because of the possibility to obtain a true non- 
equilibrium situation when two different temperatures coexist on the same manifold. On the other hand, black hole 
physics in spacetimes with positive cosmological constant appear to be of direct physical interest, because the present- 
day measurements of cosmological parameters confirm the presence of a small positive cosmological constant, which 
implies that dS backgrounds are the real black hole backgrounds to be taken into account for physical considerations. 
In this paper, we consider some relevant aspects of the physics of massive quantum Dirac particles on dS black 
hole backgrounds. We first show that, as expected, the Dirac Hamiltonian is well behaved in the sense that its 
self-adjointness can be ensured without imposing any boundary condition. We also determine, by means of spectral 
analysis, two relevant physical properties: there is no mass gap in the spectrum, even if the particles are massive, and 
there exist no quantum bound state for charged particles around a charged black hole, in contrast to classical physics. 
The latter property amounts to the absence of normalizable and time-periodic solutions of the Dirac equation on the 
background of a non-extremal Reissner-Nordstrom-dS black hole, in full agreement with the recent literature on this 
topic [3, H, 0, H, Q • Furthermore, we show that this holds true also in the extremal case, due to the prominent role 
of the cosmological event horizon, as well as in the so-called Nariai case and in the ultracold ones. 
In the second part of the paper, we also take into account the problem of pair-creation by a charged black hole. This 
is a long-standing topic in the framework of quantum effects in the field of a black hole, as old as the Hawking effect 
but still different in its origin 0, [§[ . The latter can be brought back to vacuum instability in presence of an external 
field (see e.g. 0, E3] and, in the recent literature [ll|, El, EH)- It is shown that the presence of level crossing, i.e. 
of overlap between positive continuum energy states and negative continuum energy ones, according to a criterion 
introduced by Ruffini, Damour and Deruelle [3, EE EE E3 > recently extended to include the Reissner-Nordstrom-AdS 
case [lH , is a still valid tool for investigating pair-creation of charged Dirac particles even in presence of a positive 
cosmological constant. We point out that this method is equivalent to the ones commonly exploited in order to 
investigate instability properties of the vacuum [ill EH , even if the criterion of level-crossing seems to be specific of 
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the above references iiii. 

A special attention is focused on special cases, like the Nariai and the ultracold ones, for which an estimate in WKB 
approximation of the transmission coefficient related to the process of pair-creation is provided. 

This work, together with the analogous one concerning the Dirac equation on the background of a Reissner-Nordstrom- 
AdS black hole, completes the analysis of the process of pair-creation by a charged black hole in presence of a 
cosmological constant, and in this sense it also extends the analysis on the background of a Reissner-Nordstrom 
black hole 0, 0. We recall that, in spite of the fact that dS and AdS differ for a change of sign in the cosmological 
constant, very different manifolds and very different physics occur on these backgrounds. We mention for example the 
occurrence of closed timelike curves in the AdS case, a problem which can be overcome by passing to the universal 
covering, but at the price to deal with the lack of global hyperbolicity [l9l ]. On the quantum level, self-adjointness of 
the wave operators cannot be ensured in general (see e.g. [18| for RN-AdS black holes), and boundary conditions have 
to be introduced for some cases, because of a boundary-like behavior of the AdS asymptotic region. In particular, 

for ^y^fy < 5> where ju is the Dirac particle mass and A is the (negative) cosmological constant, several boundary 

conditions can be chosen (see e.g. (20| for explicit choices of boundary conditions for the Dirac Hamiltonian on pure 
AdS), and then physics is not uniquely defined. 

In the dS cases we discuss herein, no such features arise. Moreover, black holes are characterized by a single event 
horizon in the AdS case and by two event horizons in the dS one. This fact is shown to be at the root of the fact 
that in the de Sitter case there is always level-crossing, which is in contrast not only to the AdS case but also with 
the standard RN case (A = 0). This feature is then peculiar of these solutions; notwithstanding, the actual presence 
of pair-creation is to be associated with further conditions, to be related with the actual largeness of the forbidden 
region separating positive energy states from negative energy ones. 

For completeness, we recall that charged Dirac fields in the more general Kerr-Newman-de Sitter background have 
been studied with the aim to determine their quasinormal modes in [2~D |. 



II. DIRAC HAMILTONIAN IN THE CASE r+ < r c : REISSNER-NORDSTROM-DS BLACK HOLES 

We first define the one-particle Hamiltonian for Dirac massive particles on the Reissner-Nordstrom-dS black hole 
geometry (RN-dS black hole in the following). We use natural units % = c = G = 1 and unrationalizcd electric units. 
The metric of the RN-dS black hole manifold (t E M; re (r + ,r c ); Q <E S 2 ) is 

ds 2 = -f(r)dt 2 + -^-dr 2 +r 2 dn 2 



2M Q 2 A 
r r 2 3 



M is the mass and Q is the electric charge of the black hole, and A > is the cosmological constant; let us define 

nr. 

A ' 



L — a/ 4; the equation f(r) = is assumed to admit solutions r c > r + > r_ > tq; then one obtains 



r c is the radius of the cosmological horizon, r + is the radius of the black hole event horizon, r_ is the radius of 
the Cauchy horizon. Moreover, due to the actual lack of a term proportional to r 3 , one has ro = — (r c + r + + r_ ) . 
The above reparamctcrization of the metric amounts to implementing the following relations between r c ,r + ,r_ and 
M, Q, L: 

L 2 = r c (r + + r_ + r c ) + r\ + r 2 _ + r + r_ 
2L 2 M = r 2 r + + r c r\ + 2r c r + r_ + r 2 r- + r c r 2 _ + rir_ + r+r'i 
L 2 Q 2 = r c r + r-(r c + r + + r-). 

It is not difficult to show that four real zeroes of f(r) = exist (and three are positive) if and only if the following 
conditions are implemented: 

o < Q 2 < § (3) 

M extr <M< M max , (4) 
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where 



M Rm t.r — 



3\/6 



is the mass of the extremal black hole with r_ 



1- A 1-12 



and 



Q 2 



1 - 12 



91 

' L 2 




1-12% (2-^1-12^ 



(5) 



(6) 



is the mass of the black hole with r + = r c (see sect. IIV[) . See [Hj] for the analysis of the more general Kerr-Newman-dS 



case. 



The vector potential associated with the RN-dS solution is = (— Q/r, 0, 0, 0). Spherical symmetry, as usual, 
allows to separate variables [13, [3, HE, HH and to obtain the following reduced Hamiltonian 



Hred = 



(7) 



where f(r) is the same as in (JXJ) , k = ±(j + 1/2) G Z — {0} is the angular momentum eigenvalue and \i is the mass of 
the Dirac particle. The Hilbert space in which H re( i is formally defined is the Hilbert space L 2 [(r+,r c ), l//(r) dr] 2 of 

the two-dimensional vector functions g = ' 1 



.92 



such that 



dr 

W) 



(Mr)| 2 + |<? 2 (r)| 2 )<oo. 



As a domain for the minimal operator associated with H re d we can choose the following subset of 
L 2 [(r + , r c ), 1/ f(r) dr] 2 : the set C^°(r + , r c ) 2 of the two-dimensional vector functions g whose components are smooth 
and of compact support [27[. It is useful to define a new tortoise-like variable y 



1 



and then one obtains i/gl, with y — > oo r 



where 



dr /(r) 
and y — > —oo r 
Hred = D + V(y) 



(8) 



The reduced Hamiltonian becomes 



(9) 



Do 



o a 



-cL 



and 



V(r(y)) 



e Q 



k^- 



The Hilbert space of interest for the Hamiltonian ([9]) is L 2 [R,dy} 2 . We check if the one-particle Hamiltonian is well- 
defined in the sense that no boundary conditions are required in order to obtain a self-adjoint operator. This means 
that we have to check if the reduced Hamiltonian is essentially self-adjoint; with this aim, we check if the solutions of 
the equation 



H re d g — A g 



(10) 



are square integrable in a right neighborhood of y = — oo and in a left neighborhood of y — +oo. The so called Weyl 
alternative generalized to a system of first order ordinary differential equations ([271]. theorem 5.6) can be applied, in 
particular if in a right neighborhood of y = — oo at least one solution not square integrable exists for every A S C, then 
no boundary condition is required and the so-called limit point case (LPC) is verified; if instead for every A G C all the 
solutions of (Hred — A)<7 = lie in L 2 [(— oo, c), dy] 2 , with — oo < c < oo, the so-called limit circle case (LCC) occurs 
(and boundary conditions are required). Analogously one studies the behavior of solutions in a left neighborhood of 
y = oo. The Hamiltonian operator is essentially self-adjoint if the LPC is verified both at y = —oo and at y — oo (cf. 
[27l | . theorem 5.7). In the case at hand, we can refer to corollary to theorem 6.8 (p. 99) of [27| . both for y — > — oo 
and for y — > oo. Thus, the Dirac operator defined on C Q yo (r + , r c ) 2 is essentially self-adjoint on the RN-dS black hole 
background. 
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III. QUALITATIVE SPECTRAL PROPERTIES AND TIME-PERIODIC SOLUTIONS IN THE CASE 

r+ < r c . 

We first show that the essential spectrum of the unique self-adjoint extension of the Dirac Hamiltonian [still indicated 
with H rec i] coincides with K both in the case of non-extremal black holes and in the extremal case. This feature is 
expected in presence of a black hole horizon and is well-known in the case of scalar particles [28|], and also verified 
in the case of Dirac particles on Kerr-Newman black hole manifold (see e.g. We confirm that it is verified 

also in the present cases. From a physical point of view, it also implies that there is no room for isolated eigenvalues, 
and then that there is no "standard" bound state (in the sense that a charged particle with charge opposite to the 
charge of the black hole cannot form a bound state which is analogous to the bound state an electron forms around 
an atomic nucleus). Moreover, a finer analysis allows also to conclude that, both in the non-extremal case and in 
the extremal one, the point spectrum is empty, and then no quantum bound state, i.e. no possibility to obtain a 
normalizable time-periodic solution of the Dirac equation exists. 



A. Essential spectrum 



One expects that, in presence of an event horizon, i.e. of a so-called ergosurface, the mass gap vanishes and that the 
continuous spectrum includes the whole real line. We recall that qualitative spectral methods for the Dirac equation 
(see e.g. [IS H3]) have been applied to Dirac fields on a black hole background in [!, [2§|. In order to verify this 
property, we adopt the decomposition method 27]. We split the interval (r + , r c ) at a inner point r\ and then consider 
the formal differential expression ([7]) restricted to the sub-intervals (r + ,ri] and [ri,r e ). Roughly speaking, we refer 
to the aforementioned expressions as to the "restriction of the Hamiltonian H re d to the interval (r+,n] and to the 



interval [ri, r c )" and write e.g. H, 



red\[n,r c ) 



for the latter. We limit ourselves to consider the latter restriction, which is 



relative to the novel feature of space-time, with respect to previously discussed cases, represented by the cosmological 
horizon. In the tortoise-like coordinate y one finds a potential P such that 



P 



-v7m + 



e Q 



V7m + ^ 



and it holds 



lim P(y) = P Q = 



$ c 
$ c 



which is in diagonal form and whose eigenvalues coincide. We apply theorem 16.6 p. 249 of Ref. [271 ] . which implies 
that, if v+, with < v +1 are the eigenvalues of the matrix Pq, then {R — {v- 1 C cr e (-Hr-e<j | (r-i) ,oo) ) if 



lim — 

y^oo y 



dt\P(t) - Pol = 0, 



(11) 



where | • | stays for any norm in the set of 2 x 2 matrices (we choose the Euclidean norm). In our case one has to find 
the limit as y — > oo for the following expression: 



r{y) 



dr 



1 



1 



h(r) y/r c - r 



2/j 2 h{r) 



2(<S> 2 + (r c 



1 



r) + k 2 h{r)) — 



(12) 



where we put h(r) — • -Both in the non-extremal case and in the extremal one, the above integral is finite as 
r — > r c i.e. as y — > oo, and then the limit is zero (we recall that the difference between non-extremal case and extremal 
one from this point occurs when studying the limit as r — > r+, i.e. as y — > — oo. In the extremal case r + = r_ the 
corresponding integral diverges but a trivial use of the PHospital's rule allows to find that the aforementioned limit 
is still zero). As a consequence, we can state that 



a e (H re d) — 



(13) 



A completely analogous conclusion can be stated for the restriction to (— oo,ri), and again the essential spectrum 
contribution one finds is M both in the non-extremal case and in the extremal one. 
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B. Absence of states of the point spectrum 



Qualitative spectral analysis of the reduced Hamiltonian in the non-extremal case can be implemented by means 
of theorems in [3l[ or also in [27| , In [22j a proof was given for the more general case of the Dirac equation in a 
Kerr-Newman-de Sitter black hole background, again in the case r + < r c . For the sake of completeness, we sketch 



the strategy and also provide some details involving some differences with respect to [221] . 

We note that, given a decomposition point r\ € (r + ,r c ), we can introduce the fol- 
lowing self-adjoint operators Uhor and H c on the respective domains D(Hh or ) — {g 6 
L 2 [(r + , ri), dr] 2 , g is locally absolutely continuous; g\{r\) = 0; Hhor9 € L 2 [(r + , ri), l/f(r) dr] 2 }, and 

analogously D(H C ) = {g G L 2 [{r\, r c ), l/f(r) dr} 2 , g is locally absolutely continuous; gi (ri ) = 0; H c g e 
L 2 [(n, r c ), l/f(r) dr] 2 }. According to the decomposition method applied to the absolutely continuous spec- 
trum, one has a a c{H re d) = cr ac {Hhor) U cr ac (H c ) (cf. e.g. [I2])- Theorem 16.7 in [27| allows to conclude that, in the 
non-extremal case, Hh or has absolutely continuous spectrum in R — $ + , where — ^ is the electrostatic potential at the 
black hole event horizon, and that H c has absolutely continuous spectrum in M — $ c , where — is the electrostatic 
potential at the cosmological horizon. Moreover, for eQ > it has to hold <& c < $ + , and <I> C > $ + for eQ < 0, due 
to the inequality r c > r+. In any case, $ c ^ $+ occurs, and this is an interesting fact in the light of the study 
of the pair-creation process, as we shall see in the following section. As to the spectral properties of the reduced 
Hamiltonian, one can easily infer that the spectrum is absolutely continuous in M (indeed, the above analysis allows 
to conclude that the spectrum is absolutely continuous in R — {{$ c } H {$+}} but of course the latter set coincides 
with M). 

As to the extremal case, one can again refer to theorem 16.7 in (2?| for H c and to theorem 1 in [3l[ for Hhor to 
conclude that the spectrum is absolutely continuous in K — {{^c} H {$+}}. Again, the latter set is K. 

IV. NARIAI SOLUTION 

We take into consideration the special case of the so-called charged Nariai solution [H[ , which is a black hole 
solution with r_ < r + = r c . As known, the metric U) is no more valid and a suitable transformation is necessary. It 
can be shown that the manifold can be described by 

ds 2 = I(-sin 2 (x)# 2 + d X 2 ) + 4(^ 2 +sin 2 (0)# 2 ) (14) 
A B 



with ^6lj6 (0, tt) and where B = f 1 — y 1 — 12^r ) and A = — B are constants such that < 1 [32, [3 



We note that there is no warping factor in the metric between the "radial" part and the S 2 part. For an electrically 
charged black hole we can choose Aj = — 
We study the Dirac equation as in [23|, [2^ 
Dirac matrices such that {7i,7j} = ^gif 



We study the Dirac equation as in [23|, [24| . With the same notation as in [23| , we introduce the so-called generalized 



sin(x)„ VJ . 

7o = — t=^7o 7 = — ^-^70 



A sin(x) 
7i = —j=l\ 7 1 = v^47i 



I A 

72 = ^=72 7 2 = V S7; 



2 



sin(0) , a/B , v 

73 = -7^-73 7 = ^7^:73, (15) 
VB sm(0) 

where 7,, i = 0, 1, 2, 3 are the usual Dirac matrices in Minkowski space. The Dirac equation is 

[7 fc (9 fc -r fc )- A1 ]* = o, (16) 



with 



Tfc = -77 J '(^ 7j - 7 /r^) + ieA k . (17) 
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One finds the following non- vanishing Christoffel symbols = cot(x),rj = sin(x) cos(x), T 2 3 

— sin(0) cos(0), = cot(0). Then, due to our choice for A*, we get r = — ^ cos(x)7o7i + ieA ,Ti — 0,r 2 
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0,T3 = \ cos(6')7273. Then the Dirac equation becomes 



'A /—I 



VI 72 9, 



1 



cot(e) 



■73 



sin(0) 



da 



* = 0. 



(18) 



By posing ^ — (sin.(x)) -1 ' 2 (siii(0)) _1 ' 2 £, we eliminate the terms proportional to cot(0) and to cot(x) in the previous 
equation. We now consider a static solution with £ = cxp(— iu>ip)r)(x, 0, <ft). Then a trivial manipulation of the Dirac 
equation leads to the following eigenvalue equation: 



with 



H 



H T) — LOT), 



'B 



-17071 sin(x)<9 x - eA h + \j -r sin(x)7i^ + 270-7= sin (x). 



I4 stays for the 4x4 identity matrix and K is the following operator: 

K = -i7i7o72<9<? - «7l7o73 



1 



sin(0) 



(19) 



(20) 



(21) 



which commutes with H and whose eigenvalues are [23|, [24j . By restricting H to eigenspaces of K and by 

choosing 



7o 



lJi 2 u 2 
2 -tl a 



7i 



<U 2 112 

la ©2 



(I2 is the 2x2 identity matrix, 2 is the 2x2 zero matrix), we obtain the reduced Hamiltonian 

H k = -sin(x)3 x 



J>2 — JL2 
I 2 ©2 



eQ — cos(x) 



li 2 <U* 2 

3>2 h 



= hk® I2, 



2 J12 
2 2 



/' 



sin(x) 



112 ^2 

>2 -I 2 



where 



The coordinate transformation 



eQf cos(x) - 7)= sin(x) sin(x)<9 x + y f sin(x)fc 
- sin(x)3 x + sin(x)fc eQf cos(x) + 7^ sin(x) 



.t = log(tan(— )) < — > x — 2 arctan(exp(x)) 



is such that x £ M. and, furthermore, hk becomes 



d x 
-d x 



where 



P(X) 



eQf cos(x) - sin(x) ^ f sin( X )fc 

I sin(x)fc eQf cos(x) + 7)= sin(x) 



(22) 



(23) 



(24) 



(25) 



(26) 



hk is formally self-adjoint in L 2 [E, dx} 2 and it is essentially self-adjoint in Cfi°(l&) 2 , as follows from corollary to theorem 
6.8 (p. 99) of [27| (the limit point case occurs both at x = —00 and at x = 00). It is easy to show that the essential 
spectrum of hk coincides with R and the same is true for the absolutely continuous spectrum. The latter claim can 
be checked by following the ideas displayed in sect. IIII Bl See Appendix lAl for more details. 
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V. ULTRACOLD CASE 

There is still a sub-case to be taken into account. It corresponds to the so-called ultracold case [Hj], where the 
three horizons coincide: r_ = r + = r c . Also in this case the metric (JXJ) is no more valid, and a suitable limit has to 
be considered [32| . Actually, one can introduce two different metrics for the ultracold case. As a consequence, also 
our analysis is split into two parts. 



A first metric 132] is 



ultracold I 



. x ^2 +dx 2 + ±_ {d6 2 



sin 2 (8)d(f> 2 ), 



(27) 



with x € (0, oo) and tp G R. One gets = Tq = x, r| 3 = — sin(0) cos(8), T% 3 = cot(0). The electromagnetic field 
strength is F = VAx^X A dtp, and we can choose A — -^x 2 and Aj = 0, j = 1, 2, 3. We introduce 



70 = Xlo 

71 = 7i 

72 = V2T 72 

sin(( 

73 



/2A 
sin(0) 

■73 





1 . 


7° 


= 7o 




X 


7 1 


= 7i 


7 2 


= \/2A7 




V2A 


7 3 


sin(0) 



73, 



(28) 



and then we obtain r 



■^7o7i 



■ieAo,Ti = 0,T2 = 0,T3 = | cos(0)7273. Calculations which are strictly analogous 



to the ones performed in the Nariai case (with VE" = exp(— itvip) 
lead to the following reduced Hamiltonian: 



and the variable change x = ex P( x )j & € 



ft* = 



— exp(2a;) — ^exp(x) d x + \/2Afcexp(x) 
— d x + V2Ak exp(x) exp(2x) + \i exp(cc) 



(29) 



As in the Nariai case, hk is formally sclf-adjoint in L 2 [R, dx] 2 and it is essentially self-adjoint in C^°(R) 2 , as follows from 
corollary to theorem 6.8 (p. 99) of [27J. The analysis of the spectrum can be pursued by means of the decomposition 
method applied to the absolutely continuous spectrum, and one can again conclude that the absolutely continuous 
spectrum of the self-adjoint extension of hk on M coincides with the whole real line. See Appendix IBl for details. 



B. ultracold II 

The second allowed metric [32] for the ultracold case is 

ds 2 = -dt\? + dx 2 + ^(de 2 + sin 2 (0)d0 2 ), (30) 

with x € R and ip € M. One gets T§ 3 = — sin(0) cos(0), = cot(0). The electromagnetic field strength is F = 
—yXdip A dx, and we can choose Ao = VXx and Aj = 0, j = 1, 2, 3. We introduce 



7o 


= 7o 


7° 


= -7o 


7i 


= 7i 

= V2A 72 


7 1 


= 7i 


72 


7 2 


= \/2A72 


73 


sin(0) . 
= V2A 73 


7 3 


V2A _ 
~ sin(0) 73 ' 



(31) 
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and then we obtain Tn 



ieA^T 1 = o,r 2 = o,r 3 



exp(— iuitp) 



003(^)7273. Again calculations as above (with \& 



C) lead to the following reduced Hamiltonian: 



hk 



-e^/Xx - [i d x + V2Afc 
-d x + \/2Ak -e\fkx + (i 



(32) 



Also in this case, hk is formally self-adjoint in L 2 [M.,dx] 2 and it is essentially self-adjoint in Cq°(R) 2 . As in the 
previous case, the decomposition method applied to the absolutely continuous spectrum allows to draw the conclusion 
that the spectrum of the self-adjoint extension of the Hamiltonian (|32[) is absolutely continuous and coincides with 
K. See Appendix [C] for details. 



VI. PAIR CREATION AND LEVEL-CROSSING 

We follow the Ruffini-Damour-Deruelle [HI, [TH, [H, [l7| approach, which was summarized in a previous paper [l8| . 
Herein, we limit ourselves to recall some very basic properties, focusing on the RN-dS case (the other cases can be 
dealt with analogously). In this approach one introduces effective potentials E^(r) for the positive and negative 
energy states respectively; they represent the classical turning points for the particle motion and lead to the definition 
of the so-called effective ergosphere. These potentials enter the Hamilton- Jacobi (HJ) equation for a classical particle. 



They can be interpreted also at the quantum level, as in 17J. In particular, they indicate the regions of level-crossing 



between positive and negative energy states [TEI, UH • In the case of the Dirac equation, it is known that the HJ equation 
corresponds to a WKB approximation to the Dirac equation at the lowest order [3J, [35[ ■ Variable separation in the 
quantum case allows to obtain an obvious improvement of the semi-classical formulas, amounting in replacing the 
classical value of the angular momentum with the quantum eigenvalues of the correspondin g q uantum operator [l8l |. 
We limit ourselves to recall one of the main point of our analysis of the Dirac Hamiltonian in [18[ . The key-observation 
resides in the following fact: if one consider the potential term in the Dirac Hamiltonian 



V(r) = 



Pii(r) P12M 
P2i(r) P22(r) 



and formally calculates the eigenvalues of the above matrix, which are found by solving 

(p u (r) - A)( P22 (r) - A) -pi 2 (r)p 2 i(r) = 0; (33) 
then, defining S(r) = y/ Jpvjr) + p 22 (r)) 2 - 4pn(r)p 2 2(r) + %>i 2 (r)p 2 i (r) one finds 



and moreover one gets 



\ ± (r) = ^(p 11 (r)+p 22 (r)±S(r)), (34) 



\ ± (r)=E±(r), (35) 



i.e., the semi-classical potentials introduced by Damour-Deruelle-Ruffini coincide with the eigenvalues of the matrix 
potential term in the Dirac Hamiltonian. Moreover, the square of the classical angular momentum term is replaced 
by the square of the eigenvalues fc = ±(j + 1/2) for the quantum angular momentum, and one obtains 



A. Level crossing in the RN-dS case 

Level-crossing amounts to the presence of overlap between the range of Eq and the range of Eq , signalling the 
possibility of a tunneling between positive energy states and negative energy ones; the latter phenomenon is in turn 
interpreted as pair-creation at the barrier potential, and is strictly related to the Klein paradox (incidentally, it could 
be called for this reason also Klein effect [29j). There is a peculiar property in the case of RN-dS black holes: due to 
the inequality <I> C ^ $ + , an overlap is always present, indeed one gets 

Iimf^(r)=$ +J (37) 
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and 

lira E±(r) =$ c . (38) 
r — ►rc 

Assuming for dcfinitcncss eQ > 0, one finds 

£ -(r + ) = $ + >$ c = £ + (r c ), (39) 

which proves the above claim. Moreover, it is easy to realize that level crossing occurs for energy uj such that 

<1> C = mmE^(r) < uj < max£ "(r) = <f> + . (40) 

For eQ < the overlap still exists, being E^(r c ) = $ c > $ + = £g~(r + ), and ([4*0)1 changes accordingly: $ + < cj < $ c . 
It is to be immediately pointed out that this overlap is not enough to conclude that a sensible pair-creation process 
occurs in the given background. Indeed, the potential barrier to be overcome by the negative energy particle can be 
as large as the whole external region of the spacetime. As a consequence, no efficient process can be expected on this 
ground in general, and further conditions taking into account the effective largeness of the barrier have to be looked 
for. Notice also that Eq , as a function of [i 2 , is increasing, and the same is true for its dependence on fc 2 ; both these 
properties are reversed in the case of Eq , which is in fact decreasing. As a consequence, one expects a more difficult 
level-crossing for increasing values of fi 2 or k 2 . 

We start giving sample-examples for a non-extremal black hole manifold. In figure 1, we keep fixed the geometric 
background parameters L, M, Q and the charge e of the Dirac particle, and consider two sample values of the mass. On 
the left, we find the former phenomenon described above, i.e. a level-crossing with a barrier as large as ~ r c — r + . On 
the right-hand case, which is obtained by considering a lower fermion mass, we obtain instead a level-crossing associated 
with a much smaller extent of the barrier. The latter case is expected to be involved in a effective phenomenon of pair 
creation at the barrier. In figures 2 and 3 a further non-extremal case is displayed, with ~ 1.01 (to be compared 
with ^f- ~ 10.4 of the example displayed in figure 1). In figure 3 we show the details of the potentials near the 
horizons. In figure 4, an extremal case is also shown. In figure 5 we display the so-called lukewarm case [321 13(| . which 
is such that the same temperature occurs in the case of the cosmological horizon and of the black hole event horizon, 
still with r + < r c (this happens for Q = M). It is displayed for the sake of completeness, even if from the point of 
view of the given phenomenon no peculiar behavior is expected with respect to the cases explored in figures 1 and 2. 
Notice that, from a physical point of view, the phenomenon of pair creation by a charged black hole has been related 
to the Schwinger calculation of pair creation by a homogeneous electric field [3] . The highest value of the electrostatic 
potential, which corresponds to the highest intensity of the electrostatic field, occurs near the black hole event horizon, 
hence one could naively expect that the standard condition > fi 2 , i.e. 

*)">4. (41) 
r + J e z 

which is enough for the standard Reissner-Nordstrom case @, H, [HI, is also qualitatively relevant in the present one, 
at least in some approximation. We recall that, in the case of the lightest known charged particle, i.e. the electron, 
one has — ~ 10 -21 . The above condition is not necessary for the existence of a sensible level-crossing. A proof of 

IT 2 

its sufficiency in the extremal case can be given under the hypothesis t- < 1, which allows to neglect the angular 

r + 

momentum contribution in the potentials, and also for r + <C r c . In the extremal case, it is easy to show that the 

condition to be satisfied is that -§^-{ r +)-^r( r +) > 0, which means that both potentials are increasing or decreasing 
near r = r + . One has to impose 



dE Q dE Q 1 

( r +) -r~0+) - 

dr dr 



2 



<&i-7%(i-^](i + 3^]U 



V- 



r +\ A o r +\ ( 2 , k 



',.2 



..2 



> 0, (42) 



and then, under the above hypotheses t«1 and r+ « r c , one obtains 

r + 

dE+ , s dE~ , , 1 / , , r 2 \ 

^(^)^ + )^(^-^)>0, (43) 

for which, being ^ < 1, the aforementioned condition (|41[) is sufficient. In figure 6 we show an example of extremal 
black hole where a significant level crossing occurs but condition (|41[) is not satisfied. 

Explicit evaluations of the transmission coefficient which is related to the pair-creation phenomenon (discharge) can 
be given e.g. in a WKB approximation, as pointed out in the original literature @, BEH- See also [l8[ for a short 
summary. We do not delve into quantitative evaluation herein in the RN-dS case but limit ourselves to some estimates 
in the cases which will be analyzed in the following subsections. 
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FIG. 1: On the left, we display the level-crossing in the case of a non-extremal RN-dS black hole, with L = 1000, M = 50, Q = 
30, n = 1, e = 1, k = 1. The particle and the black hole have charges with the same sign. One finds r+ ~ 90.842, r c ~ 946.214, 
r_ ~ 9.999 and ro ~ —1047.056. The upper straight line represents eQ/r+, the lower is eQ/r c . The upper potential is Eg(r), 
the lower one is Eq(t). Level-crossing occurs, but the potential barrier is as large as the whole spacetime region at hand. On 
the right, the only change with respect to the previous figure stays in the smaller value \i = 0.01 of the fermion mass. The 
upper straight line represents eQ/r+, the lower is eQ/r c . The upper potential is Eq(t), the lower one is Eq(v). Level-crossing 
occurs in this case with a much smaller extent of the potential barrier with respect to r c — r + . 





99 99.2 99.4 99.6 99.8 100 



FIG. 2: Le vel-cro ssing in the case of a non-extremal RN-dS black hole having L = V31791, M = 1193005/31791 and Q = 
330v / T90/v / 31791, which are such that r c = 100, r+ = 99 and r_ = 10. Moreover, we choose e = 1, k = 1, and the same sign 
for the black hole charge and for the particle charge. The upper straight line represents eQ/r+, the lower is eQ/r c . The figure 
on the right displays the potentials for fj, = 1, and shows that a very large potential barrier occurs. In the figure on the right 
one has /i = 0.01 and a very narrower potential barrier. Note that in the latter case eQ/r+ > fj, holds, whereas in the former 
the opposite inequality is implemented. 



B. the Nariai case 

The potentials (x) in the Nariai case are 



Set (X) = eQj cos(x) ± \J !j + k^jsmix). (44) 

Also in this case level-crossing is always present, being Eq(tt) < £^7(0) for eQ > and Eq(0) < Eq(tt) for eQ < 0. 
Level-crossing occurs for energies u> such that Eq(tt) < lu < Eq (0) in the former case and for Eq (0) < uo < Eq(tt) in 
the latter one. Again, the largeness of the potential barrier depends on the choice of the parameters. See figure 7. 

An estimation of the transmission coefficient can be given in the WKB approximation; one obtains [To] ] 



\T™ KB \ 2 = exp(-2 [ dx^), 

J barrier 



(45) 
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FIG. 3: Level-crossing in the case of a non-extremal RN-dS black hole having L = V31791, M = 1193005/31791 and Q = 
330vl90/v3l79T, as in the previous figure, with e = 1, k = 1. We display on the left the presence of a bump near the black 
hole horizon in the case of Eq . Analogously, on the right we show the behavior of the potentials very near the cosmological 
horizon. 
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FIG. 4: Level-crossing in the case of an extremal RN-dS black hole, with L = 1000, Q = 10 and then M ~ 9.999499 (see eq. 
©). Particle parameters e = 1, k = 1 are kept fixed, whereas it holds \i — 1 on the left and \i = 0.01 on the right. Level-crossing 
is more effective in the latter case, and it occurs without any bump near the black hole event horizon. 
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FIG. 5: Level-crossing in the case of a lukewarm RN-dS black hole, with L = 1000, M = 50, Q = 50, e — 1, k — 1, and with 
fi = 1 on the left and fx = 0.01 on the right. One finds r+ ~ 52.786 and r c ~ 947.214. Level-crossing is qualitatively similar to 
the one displayed in figure 1. 
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FIG. 6: Level-crossing in the case of an extremal RN-dS black hole, with L — 50, Q = 14.4336845607765725, which are such 
that M = 13.608225263871805121, r+ ~ 20.380115 and r c ~ 20.476977. With /j, = 0.01, e = 0.01, it = 1 one gets $+ ~ 0.708/i, 
which violates (|4ip . 





FIG. 7: Level-crossing in the case of a Nariai solution. The potentials E^(\) are plotted, with L = 1000, Q = 80. Particle 
parameters e = 1, k = 1 are kept fixed; particle mass is chosen to be ^ = 1 on the left and /i = 0.01 on the right. Also in the 
latter case, a bump is present for E£ near x — an d a hollow occurs for Eq near x = 7r - 

where x stays for the coordinate defined in J24])) and where we have stressed the dependence on the energy u> of 



Z -' = ( \ k " + 1") Sin2 ^^)) ~ ^ ~ eQ J cos ^))) 2 - 



(46) 



Let us introduce: 



and also 



A A 



E„ 



if. 



(47) 



(48) 



which corresponds to -4 times the maximum value for the modulus of the electrostatic field. Note that positivity of 
Z u requires to 2 < n\ + e 2 Ef n . We obtain 
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WKB\2 



exp -27r|e|J5 m (Wl 



1) , 



(49) 
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which does not depend on oj. See Appendix [Dl for more details. At the leading order as /«? <§C e 2 £'^ 1 one is also able 
to recover the approximation 

K™|*~«p(-^), (50) 

which shares a nice resemblance with the WKB estimate of the transmission coefficient related to the pair creation 
process in a uniform constant electrostatic field. Cf. e.g. [lj|. 



the ultracold cases 



We obtain 



Ef( X ) 



X~ ± \/2Ak 2 + /i 2 x 



(51) 



in the case of the first kind of ultracold solution (|27[) . and level-crossing requires that u> < for e > and u> > for 
e < 0. In the case (1301) one finds 



Efj(x) = -eVXx ± \/2Ak 2 + /x 2 , 



(52) 



It is evident that for any weR one obtains level-crossing. In figures 8 and 9 level-crossing is displayed for both the 
ultracold I and the ultracold II cases. 




FIG. 8: Level-crossing in the case of the ultracold I metric. We choose A = 0.01, k = 1, e = 1 in both cases and n — 1 on the 
left and fi = 0.01 on the right. We display only a part of the full plot (but qualitatively all relevant information is given). The 
potentials converge both to zero as x ~ * and to — oo as \ ~~ * +°o. 

As to a WKB approximation for |TL,| 2 , we get in both cases 

which again is independent from uj. Notice that, by keeping into account that for the electrostatic field one finds 
E = \/A) and with the replacement [i 2 i— > fi\ = /i 2 + 2Afc 2 , one obtains again a formula which is very similar to the 
one which is associated with the description of the pair creation in a uniform constant electric field in flat space-time 
in the same approximation, and this time no requirement about the smallness of the ratio between /i 2 and \e\E is 
imposed. 

A deeper analysis for the special cases ultracold II, ultracold I and Nariai is in progress [37l |. 

Some considerations about the problem of the choice of the quantum state playing the role of vacuum are addressed. 
If one were to assume that the positive and negative frequencies associated with the Hamiltonian define the vacuum, 
one would end up with the so-called Boulware vacuum, which is viable as the real vacuum only in the ultracold II 
case, where the background temperature is zero [36]. For a Reissner-Nordstrom-de-Sitter black hole background, a 
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FIG. 9: Level-crossing in the case of the ultracold II metric. We choose A = 0.01, k = 1, e = 1 in both cases and fi = 1 on the 
left and \x — 0.01 on the right. We display only a part of the full plot, which shows of course a linear behavior. 

further difficulty arises due to the presence of both a cosmological temperature and of a black hole temperature in 
the non-extremal and non-lukewarm cases, involving a true non-equilibrium situation. A simpler case is the extremal 
one, because of the occurrence of a single temperature for the given manifold, and the same considerations can be 
made in the case of the lukewarm solution and in the Nariai case. 

We are not aware of a rigorous construction for quantum field theory on the given backgrounds. One could expect that, 
in presence of a single non-zero temperature, suitable analyticity requirements for the fields on the extended manifold 
can lead to the thermal state as in [7[, and that "heating up" the Boulware vacuum (as it can be rigorously done in 
the case of a scalar field on a Schwarzschild black hole background [38] ) , taking into account the complication of the 
level-crossing displayed above, could be a viable solution. The instability associated with the pair-creation process 
induced by the presence of the electrostatic field generated by the black hole still remains, and gives rise to the process 
of discharge we have taken into account. Thermality of the physical state modifies such a pair-creation process but 
the transmission coefficient \T\ we have calculated for a vacuum situation still plays a role, as it is shown e.g. in Ref. 
[39l ] for the case of quantum electrodynamics in flat spacetime (see also [40] ) . One obtains that for an initial thermal 
state pair-creation is still proportional to |T| 2 with a multiplicative factor depending on the temperature. We shall 
come back to this topic in 37]. The general RN-dS case is evidently more tricky and challenging, and requires a 
non-equilibrium framework. 

VII. CONCLUSIONS 

We have shown that, on the background of a charged black hole in de-Sitter space, massive and charged Dirac 
particles are described by an Hamiltonian operator which is well-behaved both on the cosmological horizon and on 
the black hole horizon. We have also inferred that in all cases the point spectrum of the Hamiltonian is empty, and 
then there is no bound state and no normalizable time-periodic solution of the Dirac equation. The presence of two 
different horizons allows a simpler analysis even in the extremal case. Moreover, the same occurrence of two event 
horizons involving different values of the electrostatic potential is at the root of the presence in any case of level-crossing 
between positive energy states and negative energy ones. This fact per se is not enough for claiming that a sensible 
pair creation effect is present on the given manifold, due to the fact that a priori the potential barrier to be overcome 
can be very large (even large as almost the whole external manifold in the RN-dS case and in the Nariai one) and 
then the effect is expected to be very suppressed. Nevertheless, in all cases examples can be found where the barrier 
is of much more reduced extent, in such a way to allow a physical ground to the pair-creation phenomenon. Some 
estimates in WKB approximation have been given for the transmission coefficient which is related to the pair-creation 
process 0, H [r| in the case of the Nariai geometry and in the ultracold ones. 
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APPENDIX A: ABSOLUTELY CONTINUOUS SPECTRUM IN THE NARIAI CASE 

We introduce a decomposition point d £ R and also the following self-adjoint operators i?_oo and ffoo on the 
respective domains Z)(iJ_ 00 ) = {g e L 2 [{— oo, d], dx] 2 , g is locally absolutely continuous; g\ (d) = 0; H- X g S 
L 2 [(— oo, d], dx] 2 }, and analogously D(H X ) = {g S L 2 [[d, oo), dx] 2 , g is locally absolutely continuous; <?i(d) = 
0; i?ooS S L 2 [[d, oo), dx] 2 }. We define P_ := lim^-oo P(x(x)) and P + :— lim. I ;^_ 00 P(x(x)), where the P(x(x)) is 
the potential (f2"6")) . and write 

P = P T + (P-P T ). (Al) 

The first term on the right hand side of eq. (|A1[) is obviously of bounded variation, whereas the latter term is such 
that |P — P+| G L 1 [[d, oo), dx] and |P — P_| € L 1 [(— oo, d], dx] respectively. Moreover, notice that 



±eQf 
±eQ| 



(A2) 



As a consequence, in both cases the hypotheses of Theorem 16.7 in [27J are implemented, and one is allowed to 
conclude that P_oo has absolutely continuous spectrum in R — {eQ^}, and that has absolutely continuous 
spectrum in R — {— eQ^}. Then the absolutely continuous spectrum of the self-adjoint extension of the Hamiltonian 
operator (j2"5|) is R. 



APPENDIX B: ABSOLUTELY CONTINUOUS SPECTRUM IN THE ULTRACOLD I CASE 

Let us introduce a self-adjoint extension -ff_oo of the formal differential expression (f2T))) on the interval (— oo,0] (0 
is the decomposition point). Notice that the potential term in (|29|) is 

p, x \ _ -^-exp(2x) - ^exp(x) V2Akexp(x) ^ 

V2Afcexp(x) ex P( 2x ) + ^exp(x) 

and it is such that lim^^-oo P(x) = O. Moreover, it is easy to show that \P{x) | <G L 1 [(— oo, 0] , da;] . As a consequence, 
Theorem 16.7 of [271 ] can be applied and the given self-adjoint extension has absolutely continuous spectrum in 
R — {0}. It is also true that is not an eigenvalue for -ff_oo, because no normalizable solution exists as a consequence 
of Levinson theorem (whose applicability is related to the property that each entry in P(x) is integrable near x — — oo; 
cf. [4l[, p. 8). Cf. also [ij for the Kerr-Newman case. Thus a a c{H-oo) = R- As a consequence (cf. e.g. 22]), also the 
absolutely continuous spectrum of the self-adjoint extension of hk on R coincides with the whole real line. 



APPENDIX C: ABSOLUTELY CONTINUOUS SPECTRUM IN THE ULTRACOLD II CASE 



Let us notice that the equation h^g = Xg, by putting g{x) 
equivalent to the following equation: 

d 



W\ (x)iii (a;) 

i / \ ), where wAx) = exp(— 2\/2Akx), is 



dx 



u(x) 



(— eVXx — a)w\(x) 

(eVAx + P)^ 



u(a;), 



(CI) 



where a = A — /i and (3 = A + a . Then, by restricting our attention to the interval [c, oo), with c > 0, 
and by applying Theorem 2 in [3l[ the result follows. Let us define, according to the notations of [3l| . 
Pi = (— eVAx + fi) exp(— 2v / 2A7cx) —: pu + pi 2 , p2 = (—eyAx — ^) exp(2\/2A£;x) =: P21 + P22, «i = exp(— 2v / 2Afcx) 
and a 2 — exp(2\/2Afcx); p\\ := — e-^/Aa; exp(— 2V2AA;a;) and p2\ '■= — e-s/Ax exp(2v / 2A/ca;). One obtains 



P11P21 



= e 2 Ax 



> and / c + °° dxy/pnp2i = +00. Moreover both J^°° dxaiJ 1 ^ and J c + °° dxa 2 \/p^ diverge 



+00 



(it is sufficient that one of them diverges [31J). Moreover, if 77 = (^j) 1 / 4 , then A := 



dx r) v /pnP2l 



^l^M i s sucn that 

|e|x 



lima;- 



A = and f± e L 1 [[c, +00), dx]. Moreover, 2±. 



a2_ £12 P22 are i ori g_ ran g e ( m th e sense that they vanish 



dx LL ' /> J' ' pu ' P21 ' pu ' P21 

as x — > +00 and their derivative is in L 1 [[c, +00), dx]). Then the hypotheses of Theorem 2 in [31| are implemented, 
which means that the absolutely continuous spectrum is R. 



16 



APPENDIX D: EVALUATION OF THE NARIAI TRANSMISSION INTEGRAL 

We need to evaluate 

dx (Dl) 



/ 

J barri 



where Z^ is given in (|46|) and x( x ) = Sarctane 31 (cf. (|24[) ). Using 

cosx(a^) = — tanhx, 
sinx(ir) = 



cosh: 



we can rewrite 



dx = I \ I ( —rk 2 + !—r) — ( djcosha; + eQ— sinhx 



barrier J barrier y / V ^ / COSh X 



If we define fi\ as in (|47|) and cj± = cj ± eQ-j, and change variable to z = e 2a; we get 



1 / ixP" 1 / 1 \ 1 u 2 

Z W = I d Z — ^--- 2 + (fi - J - - -±, (D3) 

where < z_ < z + are the two real solutions of 

uPl ( 2 1 \ or? 2 
— 4 - + I (J-k - 2 W + W - J z ~ — z = °- 

Note that such solutions exist if the discriminant of the polynomial is positive. This gives 

B 2 



The integral is indeed well defined for —y^t 2 + e 2 Q 2 ^j < uj < y fJ-f. + e 2 Q 2 ^-. Note however that, only the region 
of level crossing —\eQB/A\ < uj < \eQB/A\, corresponding to ui + oj- < 0, is relevant for computing the transmission 
coefficient. Let us first distinguish the "generic case" ut+oj- ^ from the "particular case" ui+ui- = 0. 
In order to compute the integral in the generic case, the residue method is used. Let us cut the complex plane C 
along the segment [z-,z+] C K+ (note that z_ > in this case), so defining a Riemann sheet for the square root 

f(z) = yj (z + — z)(z — z-)uj+/4. In particular we choose the phase of (z + — z)(z — Z-)u>\/& to be (modulo 4-7r) 

along the lower border of the cut, so that it will be 27r (modulo An) along the upper border. Thus, it makes sense to 
take the closed path r — Sj U c + U Sj U c_ as in figure flUl 



FIG. 10: Circuit Y around the cut in the Riemann sheet of / 

When the radii of the two circles c± are set to zero, V approaches the cut without crossing any singularity so that 
the value of the integral 

I=fdz 1 f{ - Z) (D4) 
J r 1 + z z 
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does not change. In this limit, the contributions from the circles vanish, whereas 



T 

and being the phase of / equal to on the lower cut and to it on the upper cut, we see that 

Z/j dx 



barrier 



1 + z z 2 



(D5) 



To compute the integral /, let us blow up T, without crossing the poles of the integrand (which are z = 0, — 1, oo) like 
in figure QT] 



-1* o 




FIG. 11: Blow up of the path F 

We see that (cf. figure (fTTjl) 

/c e (0) l + z z Jc4-l) l + z z Jc R {0) l + z z 

where we used c r (z) to indicate the counterclockwise oriented circle with center z and radius r. Note that if we take 
the change of variable z — 1/i in the last integral, we have [42j cr(0) — ► — c 1 / i? (0) so that 



c e (0) 1 + z z Jc„(-1) 1 + z z Jc 1/B (0) 1 + t 



-2tt* (/(0)-/(-l)-lim(z/(l/z))) 



2- [ + \/-(Ml + > + - --) 2 ) + \/- Z f I • (D7) 



To compute the square roots we need to specify their phases. This is easily done looking at the Riemann sheet. 
Indeed, if we look at the real axis, on z < z_ the phase of / 2 is 3ir (modulo 4tt) so that f(x) — — i\f(x)\ for x = 0, — 1. 
For the last root, we note that 

limz/(l/z) = lim f(z)/z 

and because the monodromy of z = oo is trivial (the phase of f(z) does not change (modulo 2n) if \z\ is very large 
and arg(z) varies by a period), z = oo is indeed a pole (and not a branch point) and this limit does not depend on 
the phase of z. Thus, we can compute it along the positive real axis. But there, the phase of f(z) 2 , for z > z + is 7r, 
so that f(z) = i\f(z)\, and finally we have 
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= 2tt 



\j^l + \^+-^-) 2 - \(\<*>+\ + |W-|) 



(D8) 

w_ I which does not 



Note that in the physically interesting case, that is when lo + lu- < 0, we have \lu + \ + = |u; 
depend on u, and reproduces exactly (j4l?]) . 

It remains only to check the particular cases, which are however easily obtained by direct integration. For example, 
for w_ = 0, so that u)+ = 2eQB/A we have 



7/2 = 



Zrj dx — 



1 

2 Jo 



dz 



1 + z V z 



2 2 



Introducing the new integration variable s 2 = — oj"j_ we get 

1 1 



7/2 = f 
Jo 



V 2 + cj 2 + s 2 



s 2 ds = 



4/i 2 +^ 2 



4^ 2 + LO Z + + s 2 



ds 



Mfc + Jw 2 



For cj + = (w_ = — 2eQ£?/A) we have 

7/2 = 



1 / 2 _ut 
Z (l + z)V^ Z 4 



dt 

Y+t V f 



2 2 



(D9) 



where we used the change of variable t = 1/z, obtaining the same integral as in (|D9[) with w_ in place of giving 
thus the same result, being \ui±\ = 2\eQB/A\. 
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